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Abstract 

In this work the dynamics of self-gravitating systems composed by dark and baryonic matter is 
analyzed. Searching for a description of this dynamics, a system of collisionless Boltzmann equations 
for the two constituents and the Poisson equation for the gravitational field are employed. Through 
the solution of these equations the collapse criterion is determined from a dispersion relation. The 
collapse occurs in an unstable region where the solutions grow exponentially with time. It is shown 
that the unstable region becomes larger if the dispersion velocity of dark matter becomes larger than 
the one of the baryonic matter. The results obtained are also compared with the case where only the 
dark matter is present. The model of the present work has a higher limit of instability and therefore, 
exhibited an advantage in the structure formation. 


1 Introduction 

The problem of structure formation has been investigated for a long time, since the pioneer work of Jeans 
[1], There are many models which explain the conditions for the beginning of the structure formation 
and when they occurred. The Jeans model describes the gravitational instability of a self-gravitating gas 
cloud and search for the conditions that small perturbations of a gas cloud grow exponentially, leading to 
the collapse of the cloud (see e.g. [2, 3, 4, 5, 6, 7, 8, 9, 10]). The simplest way to understand these criteria 
is to think in force balance arguments. The collapse occurs whenever the inwards directed gravitational 
force is larger than the outwards directed internal pressure of the gas. The parameters which are used to 
quantify this instability are the Jeans length and Jeans mass that will be discuss throughout this paper. 

In this work the dark matter [11] is present together with baryonic matter. The accepted model of 
this unknown component is the model of cold dark matter. Such model postulates the dark matter as 
weakly interacting massive particles and whose typical velocity is much lower than that of the light. As 
a consequence, the cold dark matter is an auxiliary matter in a process of structure formation. In this 
scenario, the cold dark matter leads to a hierarchical formation process where small structures are formed 
first and massive structures later. Models with hot dark matter, on the other hand, show difficulties to 
explain the galaxies formation and other structures on a small scale. The higher velocities do not allow the 
agglomeration in the necessary scale. Thus, the model which will be analyzed in this work is composed 
by a mixture of two constituents, namely, a cold dark matter and a baryonic matter subjected to a 
gravitational field. 
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For the search of the collapse criteria for our model, we invoke the system of Boltzmann equations 
for the constituents and the Poisson equation for the gravitational field. This set of equations leads to a 
dispersion relation which implies in a collapse criterion for infinity homogeneous fluid and stellar systems. 

The aim of this work is: to obtain the Jeans length and the Jeans mass, to compare the results 
with the case when only one constituent is present and to discuss the dispersion velocities ratio between 
baryonic and dark matter. 

Some works with dark and baryonic matter were analyzed in the references [3, 4, 6, 10] by taking into 
account the macroscopic balance equations, but to the best of our knowledge the use of the system of 
Boltzmann equations is new. 

Another point of view without the inclusion of a dark matter constituent is to analyze the Jeans 
instability by using a f(R)~ theory, which was the subject of investigation of the works [12, 13, 14]. 

This paper is organized as follows. In section 2, we present the set of equations, the dispersion relation 
and we explore the collapse limits through the Jeans mass and Jeans length for a system of collisionless 
baryonic and dark matter. The conclusions of the work are stated in section 3. 


2 Systems composed of collisionless dark and baryonic matter 


We start with the three equations that describe systems composed by dark and baryonic matter subjected 
to a gravitational field <f>. Here we use the indices d and b for the dark and baryonic matter, respectively. 
The evolution equations of the distributions functions of baryonic f b = /(x, v b ,i) and dark matter 
fd = f(x,Vd,t) are the collisionless Boltzmann equations (see Appendix) 
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The gravitational field must fulfill the Poisson equation: 


V“$ = 4t tG 


f b dv b + / f d dv d ) = 4t iG{p b + p d ), 


(1) 

( 2 ) 

(3) 


where p b and pd are the mass densities of the baryonic and dark matter, respectively. 

We consider that the equilibrium of the self-gravitating system is described by homogeneous and 
time-independent distribution functions /°(v b ), /°(v d ) and a potential which is only a function of the 
space coordinates, i.e., 4 > o(x). The equilibrium state is then subjected to small perturbations represented 
by plane waves of frequency w and wavenumber vector k, namely 


/(x,v b ,t) = / b °(v&) + f b exp[i(k-x- wt)],, 

/(x, v d , t) = fd (v d ) + T d exp[i(k • x - ujt )},, 
$(x, t) = 4>o(x) + 4>i exp[*(k • x - ut)], 


(4) 

(5) 

( 6 ) 


where the overbarred quantities denote small amplitudes. 

The equilibrium for a homogeneous system is achieved by Jeans ’’swindle” that allows us to make 
•Fo = 0 without loss of consistency. Hence, if we insert (4) - (6) into the system of equations of (1) - (3) 
and linearize the resulting equations, we get 


-iufb + v fe • (* k /h) - (*k<Fi) • = 0, 


(7) 
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( 8 ) 


-iufl + v d ■ - (* k$ i) ' 

-k 2 ¥i = 4 nG (Jfldv b + Jfddvd') . (9) 


Here the modulus of the wavenumber vector was denoted by k = |k|. 

By eliminating the overbarred quantities from the system of equations (7) - (9) it follows the dispersion 
relation 


k 2 + 47 tG k • 


dfS dv b | fdf° d dv d V Q 
dv b v b • k - lj J dv d v d • k - w / 


( 10 ) 


From now on we assume that the distribution functions and f d are the Maxwellians: 


fb(y b ) 


Pb -vll2al 
{2ncrl)i 


fd(Vd) 


Pd -vl/lal 
{2-Kcj 2 d )i 


( 11 ) 


where p®, p d are constant mass densities of baryons and dark matter, respectively, and a b and ad their 
dispersion velocities. 

Without loss of generality we can choose k = (k, 0,0) so that the dispersion relation (10) together 
with (11) can be integrated with respect to the velocity components v y and v z , yielding 


k 2 



Ph_ r°° x 2 e x 
a b Jo x 2 — w 2 / (2cr 2 k 2 ) 


Pd f°° v 2e v 
°dJ0 V 2 — u,2 /(2c 2 k 2 ) 


Above, we have introduced new integration variables 


dhx _ b'dx 

V2a b y \[2od 


( 12 ) 


(13) 


Unstable solutions are such that Sft(w) = 0 and w/ = Q'(w) > 0, since in this case the solutions grow 
exponentially with time. When u = iui the integrals on the right-hand side of (12) can be evaluated by 
using the following relationship (see [15] eq. 3.466) 



x 2 e 

x 2 +/3 2 



^-e^ 2 erfc(/3p), 


(14) 


where erfc is the complementary error function. By taking into account (14) the dispersion relation (12) 
with ui = iu>i reduces to 
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(15) 


In the above equation we have introduced the dimensionless wavenumber and the dimensionless frequency 
defined by 
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Note that we have used the mass density pd and the dispersion velocity ad of the dark matter in order to 
construct the dimensionless wavenumber and frequency, since the dark matter plays an important role in 
the cosmic structure formation. 

From (15) we may determine the instability regions in the plane (w*,k*), where disturbances grow 
with time. For that end we have to specify: (a) the mass densities ratio p d /p® and (b) the dispersion 
velocities ratio ad/ab . The mass densities ratio can be taken as the ratio of the densities parameters 
fld/flh today which is approximately 5.5 (see the section 22 of [16]), i.e., p° d /p^ ss 5.5. This ratio has not 
modified too much during the evolution of the Universe. For the dispersion velocities ratio we take the 
values given in the work [17] on the Milky-Way-like galaxy simulations including both dark matter and 
baryons, namely, ad/ab = 170/93 ~ 1.83. 

By considering w* = 0 in the dispersion relation (15), we obtain that the dimensionless wavenumber 
is equal to k* = 1.2694. This value can be interpreted as the ratio of the Jeans wavenumbers one referring 
to the dark matter-baryons system k db and the other only to dark matter kj = \JAnGp^/ad- Hence, the 
ratio of the Jeans wavelengths A db = 2ir/k db and A d = 27r/kj is X db /X d = 0.7877, showing that the Jeans 
wavelength of the system dark matter-baryons is smaller than the one with only the dark matter. This 
reflects also in the value of the Jeans mass, which is defined as the mass contained in a sphere of diameter 
equal to the Jeans wavelength. The ratio of the Jeans masses of the systems dark matter-baryons M db 
and only dark matter M d is 


Mf , P°d + P°b (xf \ 3 

M d j p° d V X/) 


Pd + Pb jf P>1 
Pd \l\P>b+P> 2 d 


0.5791. 


(17) 


This shows that in the model where baryonic and dark matter are present, the structures began to form 
earlier, at the time that the dark component dominated, which reinforces the fact that a smaller Jeans 
mass was required to initiate the collapse. 



Figure 1: Dimensionless wavenumber versus dimensionless frequency. Models with baryonic and dark 
matter with: (a) ad/ab = 1-83 thick line; (b) ad/ab = 1-20 dotted line; (c) ad/ab = 2.2 dot-dashed line. 
Model with only dark matter dashed line. 
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We may ask about the role of the dispersion velocity ratio on the structure formation. This can be 
seen from Fig. 1 where three different values for <Jd/&b were considered. The regions below the lines 
refer to the positive values of the contour plot of (15), while the regions above correspond to the negative 
values. The arrow in this figure indicates the direction of the increase of the dispersion velocities ratio. 
Here we may infer that by increasing the dispersion velocity ratio the unstable region also increases. As 
a consequence the ratios of the Jeans masses become: (a) Mj b /Mj = 0.8338 for <7dl<Jb — 1-20 and (b) 
Mj b /Mj = 0.4585 for Odfob = 2.2. Hence, the increase of the dispersion velocity ratio implies that a 
smaller mass is needed to begin the collapse. 

3 Conclusions 

In this work, we analyzed the dynamics and the collapse of collisionless self-gravitating system composed 
by dark and baryonic matter. This system is described by two Boltzmann equations, one for each 
constituent, and the Poisson equation for the gravitational field. The dispersion relation for this model 
was found by linearizing the equations about the equilibrium of the self-gravitating system, which is 
described by homogeneous and time-independent Maxwellian distribution functions and a potential which 
obeys Jeans ’’swindle”. The dispersion relation allows to obtain the collapse criterion of interstellar gas 
clouds composed by the mentioned constituents, resulting in the formation of structures. 

The Jeans mass always has lower values in the models where the baryonic and dark matter are present. 
This phenomena occurs because the dark matter agglomerates more easily and the baryonic matter passes 
to aggregate due to the attraction of the gravitational field generated from this initial agglomerate. This 
result demonstrates that the structures began to form earlier when the dark matter dominated. We 
also conclude that when the dispersion velocity of the baryonic particles is relatively smaller than the 
one of the dark matter, they easily aggregate, since they hardly overcome the escape velocity of a given 
gravitational field. This behavior can be observed through the Fig. 1, when dd/ub increases, the unstable 
region - where the solutions grow exponentially - increases also. 

Although the model described here is a function of two ratios - namely, of the densities and of the 
dispersion velocities of baryonic and dark matter - and the values of these ratios must be specified in order 
to analyze Jeans instability, we think that it is appropriate to analyze collisionless self-interacting systems 
of dark and baryonic matter. Here a kinetic model with two Boltzmann equations one for each constituent 
coupled with the Poisson equation was considered. It is the kinetic counterpart of the phenomenological 
model which considers the balance equations of mass and momentum for each constituent coupled with 
the Poisson equation (see e.g. [3, 6]). 

To sum up, this model with baryonic and dark matter proved to have a higher limit of instability and 
therefore, exhibited an advantage in the structure formation. 
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A The Boltzmann equation 

A state of a non-relativistic gas mixture of r constituents is characterized by the set of one-particle 
distribution functions f a = /(x, v Q ,f), (a = 1,... , r) in the phase space spanned by (x, v a ), such that 
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f(x,v a ,t)dxdv a gives the total mass of particles of constituent a in the volume element dx about the 
position x and with velocity in the range dv 0 about v a . The space-time evolution of the one particle 
distribution function f a in the phase space is governed by the Boltzmann equation, which in the presence 
of a gravitational potential $ and in the absence of collisions between the particles, reads (see e.g. [18]) 


dL 

dt 


"ja. dfa 8 f a 

^T + v a -—-V4> • —— = 0, 0 = 1, ...r, 

< 7 X a 


(18) 


where the force which acts on a particle of constituent a is only of gravitational nature F = —V$. 

The above equation follows also from a relativistic Boltzmann equation in the presence of gravitational 
fields, where a one-particle distribution function f a = /(x, p a ,t) of constituent a in the phase space 
spanned by (x,p 0 ) satisfies the Boltzmann equation (see e.g. [19]) 


i & fa 
1 dx^ 


-T! 




/V — = 0 

r'a o '- , 5 

O'Pa 


a = 1,... r. 


(19) 


Here the mass-shell condition p^Pay = rn 2 c 2 - where m a is the particle rest mass of constituent a - 
was taken into account. In the non-relativistic Newtonian limiting case p° a — > m a c , p a — > m a v a and the 
Christoffel symbol Tq 0 —> V*$/c 2 so that (18) follows. 

The equilibrium distribution function is a Maxwellian distribution of the velocities 


/°(x,v Q ,t) 


P a r -vll 

( 2^)1 


( 20 ) 


where a a = y/kT a /m a - with k denoting the Boltzmann constant and T a the temperature of the con¬ 
stituent a - is the dispersion velocity and p a the mass density of constituent a, which is defined in terms 
of the one-particle distribution function by 


Pa, = 



( 21 ) 
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